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NOMENCLATURE 

2 
A^  -  Total  body  surface  area  (DuBois)  ft 

2 
A  -  Convective  area,  ft 

c  ' 

2 
A  -  Surface  area  of  the  environment,  ft 

A  -  Radiative  surface  area  of  the  body,  ft* 

C  -  Rate  of  convective  heat  transfer,  Btu/hr 

C  -  Modified  coefficient  of  convective  heat  transfer,  Btu/hr-f t2-F1,25 
c 

CF  -  Conversion  factor  1.90633  Btu   cmHg  var 

Kcal  psia  ft2 

C  -  Specific  heat  at  constant  pressure,  Btu/lb-F 

E  -  Rate  of  evaporative  heat  loss,  Btu/hr 

F   -  Shape  factor  from  body  n  to  body  m 

g  -  Acceleration  of  gravity,  ft/hr2 

Gr  -  Grashof  number 

H  -  Height  of  an  individual,  in 

hc  -  Coefficient  of  convective  heat  transfer,  Btu/hr-f t  -F 

hfg  -  Heat  of  vaporization,  Btu/lb 

K  -  Rate  of  radiative  heat  transfer,  Btu/hr 

k  -  Thermal  conductivity,  Btu/hr-f t-F 

L  -  Characteristic  geometrical  dimension,  ft 

M  -  Rate  of  metabolic  heat  produced  within  the  body,  Btu/hr 

m  -  Mass  of  the  body,  lbm 

Nu  -  Nusselt  number 

Pa  -  Saturation  vapor  pressure  at  ambient  temperature,  psia 

Pr  -  Prandtl  number 

Ps  -  Saturation  vapor  pressure  at  skin  temperature,  psia 

R  -  Rate  of  radiative  heat  transfer,  Btu/hr 
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Re  -  Reynolds  number 

RH  -  Relative  Humidity,  % 

S  -  Rate  of  heat  storage,  Btu/hr 

Ta  -  Dry  bulb  temperature  of  the  ambient  fluid,  F 

T. &  MRT  -  Mean  radiant  temperature,  F 

MRT  *  ' 

Tg  -  Skin  temperature,  F 

AT  -  Temperature  difference  between  the  skin  and  the  ambient  fluid,  F 

dT/dt  -  Time  rate  of  change  of  temperature,  F/hr 

V  -  Velocity  of  fluid,  ft/hr 

VCO2  -  Volume  of  COo  produced,  liters /hr 

VO2  -  Volume  of  O2  consumed,  liters /hr 

W  -  Rate  of  work  output,  Btu/hr 

W  -  Weight  of  an  individual,  lb 

W  -  Fraction 

Wf  -  Subject  weight  after  1  hour,  lb 

W.  -  Subject  initial  weight,  lb 

3  -  Coefficient  of  volumetric  expansion,  F 

e-,  -  Emissivity  of  the  human  skin 

E2  -  Emissivity  of  the  surroundings 

2 

u  -  Coefficient  of  heat  transfer  by  evaporation  Btu/'ft  -hr-psia 

u  -  Viscosity  lb/hr-ft 

p  -  Density  lb/ft3 

-8         2  A 
a  -   Stefan  Boltzman  constant  0.1713  '10   Btu/hr-ft  -RH 

<|>  -  Functional  relationship  between  variables 
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INTRODUCTION 

Deep  sea  diving  and  the  prospect  of  living  and  working  on  the 
ocean  floor  (continental  shelf)  pose  several  problems  not  encountered 
in  normal  surface  living.   One  of  the  major  problems  is  that  the  com- 
position of  normal  atmospheric  air  is  no  longer  suitable  at  the  pressures 
encountered  in  the  deep.   The  oxygen  partial  pressure  (P0„)  must  remain 
within  the  limits  of  1.0  -  0.17  atmospheres.   If  it  drops  lower  than 
this,  the  symptoms  of  anoxia  become  evident  and  eventually  death  will 
result.   If  the  PO  becomes  considerably  higher,  oxygen  poisoning 
becomes  a  problem   (1) .   The  nitrogen  which  makes  up  practically  all 
of  the  remaining  air  becomes  narcotic  at  partial  pressures  above  29  psi, 
which  for  normal  air  is  a  compression  to  a  depth  of  50  ft  (2). 

Helium  lacks  the  anesthetic  effects  of  nitrogen.   For  this  reason, 
the  Navy  uses  Helium  -  Oxygen  mixtures  rather  than  air  for  its  deep 
diving  and  such  deep  sea  explorers  as  Link  of  Ocean  Systems,  Costeau 
of  Conshelf  and  Bond  of  Sealab  all  use  helium  as  their  diluent  gas. 

The  question  now  arises;  since  helium  has  different  physical  and 
thermal  properties  than  nitrogen,  how  are  the  basic  heat  transfer  and 
comfort  equations  for  the  human  body  altered  due  to  the  use  of  helium? 

The  objectives  of  this  study  were  (1)  to  determine  environmental 
conditions  which  produced  subjective  comfort  in  helium  at  one  atmos- 
phere and  to  determine  the  heat  losses  which  occur  at  these  conditions 
for  a  nude  human  subject  seated  at  rest  and  (2)  to  formulate  the 
appropriate  theoretical  equations  to  predict  the  environmental  con- 


ditions  for  thermal  balance  in  helium  at  twenty  atmospheres  pressure 
which  corresponds  to  a  depth  of  645  ft,  i.e.  ocean  exploration  on  the 
continental  shelf. 


Literature  Review 

A  review  of  published  literature  revealed  only  two  publications 
(3,  4)  in  which  a  heat  balance  study  had  been  done  on  a  human  subject 
in  a  helium  environment.  Epperson  et  al.  (3)  studied  four  male  sub- 
jects at  rest  and  during  exercise  at  atmospheric  and  reduced  pressure 
while  exposed  to  a  helium-oxygen  environment,  and  compared  the  results 
with  those  obtained  for  the  same  test  in  an  air  environment.   A  lowering 
of  skin  temperature  in  the  helium-oxygen  environment  as  compared  to 
the  air  environment  was  reported.   The  heat  loss  by  convection  was 
greater  in  the  helium-oxygen  environment  and  attempts  to  explain  the 
difference  dealt  with  forced  convection  parameters  only,  and  not  free 
convection  parameters. 

Fox  et  al.  (4)  also  studied  the  thermal  response  of  several  sub- 
jects exposed  to  an  atmosphere  of  helium-oxygen  while  resting  and 
during  exercise.   These  studies  were  concerned  mainly  with  the  differ- 
ences in  skin  and  rectal  temperatures  of  the  subjects  in  the  two  media 
as  affected  by  environmental  temperature  and  relative  humidity.   It 
was  concluded  that  the  skin  temperature  is  lower  in  He-CL  than  in  air 
and  that  the  difference  is  a  linear  function  of  the  temperature 
difference  between  the  skin  and  gas.   It  was  concluded  also  that  the 
evaporation  of  moisture  from  the  body  was  not  significantly  different 
in  the  two  media.   There  was  no  interaction  between  the  relative 
humidity,  the  thermal  properties  of  the  gas  media  and  the  evaporative 
weight  loss  rate.  The  attempt  to  explain  differences  in  convective 
heat  loss  also  used  forced  convection  parameters  only. 


THEORETICAL  ANALYSIS 

Energy  Balance 

When  man  is  in  thermal  balance  with  his  environment,  his  energy 
input  equals  his  energy  output.   This  is  true  whether  he  is  lying 
still,  walking  slowly  or  running  rapidly.   If  someone  is  subjectively 
"comfortable",  they  are,  or  should  be,  in  thermal  equilibrium  with 
their  environment. 

In  the  zone  of  thermal  equilibrium  a  person  requires  no  effort 
on  his  part  to  maintain  equilibrium,  other  than  the  normal  physiologi- 
cal mechanisms  of  vaso-motor  regulation  and  sweating.   These  mechanisms 
are  continuously  in  operation  to  maintain  the  body  at  the  desired  tem- 
perature.  If  the  body  starts  to  cool  the  blood  vessels  adjacent  to 
the  skin  constrict,  reducing  the  flow  of  blood  to  the  extremities  and 
outer  skin  layers  and  in  this  way  restrict  the  heat  flow  to  the  outer 
surfaces.   This  is  the  mechanism  of  vaso-motor  regulation  against  cold. 
If  the  body  becomes  over-heated,  the  blood  vessels  near  the  body  sur- 
face dilate  in  an  attempt  to  lose  heat  more  rapidly.   This  is  the 
mechanism  of  vaso-motor  regulation  against  heat.   A  point  is  reached 
in  the  heating  or  cooling  process  when  the  normal  vaso-motor  mechanisms 
are  not  capable  of  maintaining  the  body  temperature.  When  this  occurs 
during  cooling,  the  mechanism  of  shivering  is  employed  by  the  body  to 
generate  energy  internally  and  thus  offset  the  high  heat  loss  at  the 
body  surface.   When  the  body  begins  to  overheat,  the  sweat  glands  put 
additional  moisture  on  the  body  surface.   This  moisture  evaporates 
taking  the  necessary  heat  of  vaporization  from  the  body  and  cooling  the 
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body  in  the  process.   This  high  sweat  rate  may  be  called  "sensible 
sweating"  to  distinguish  it  from  the  low  moisture  loss  rate  which  is 
occurring  continuously.   Once  the  entire  body  surface  area  is  covered 
with  moisture,  a  maximum  evaporative  heat  loss  will  occur.   If  this 
heat  loss  is  not  enough  to  compensate  for  the  heat  gained  and/or  pro- 
duced, the  body  temperature  will  rise. 

Since  shivering  is  actually  muscular  work  which  requires  internal 
energy,  it  will  be  considered  outside  the  category  of  continuous  normal 
physiological  processes.   This  places  the  zone  of  thermal  equilibrium 
between  the  onset  of  shivering  and  point  of  maximum  cooling  by  evapora- 
tion of  sweat.   The  comfort  zone  must  lie  between  these  points,  depending 
upon  the  individual.   For  general  purposes,  maximum  subjective  comfort 
for  a  sedentary  individual  will  occur  between  a  point  just  above  the 
onset  of  shivering  and  a  point  just  below  the  onset  of  visible  sweating. 

The  zone  of  thermal  equilibrium  and  the  zone  of  subjective  com- 
fort within  it  occur  when  the  energy  gained  or  produced  equals  the 
energy  lost.   The  physical  mechanisms  by  which  this  energy  is  gained 
or  lost  can  be  represented  by  a  mathematical  equation.   This  report  is 
concerned  with  the  effect  of  environmental  conditions  on  that  equation. 

In  the  following  sections,  the  heat  balance  equation  will  be 
first  stated  in  general  terms,  and  then  analyzed  term  by  term,  as  it 
applies  to  the  experiment  of  man  in  a  helium  environment. 

The  general  heat  loss  equation  may  be  stated  as: 

M-W  =-E±K±R+C±S  |l| 

where 

M  =  rate  of  metabolic  heat  produced  within  the  body,  Btu/hour 


W  =  rate  of  mechanical  work  output,  Btu/hr 

E  =  rate  of  evaporative  heat  loss,  Btu/hr 

K  =  rate  of  conductive  heat  loss,  Btu/hr 

F  =  rate  of  radiation  heat  loss,  Btu/hr 

C  =  rate  of  convective  heat  loss  or  gain,  Btu/hr 

S  =  rate  of  storage  of  energy  in  the  body,  Btu/hr 
Metabolism  -  Heat  Production 

Energy  is  produced  within  the  body  by  the  oxidation  of  foodstuffs. 
The  rate  at  which  the  oxidation,  and  with  it  production  of  energy  takes 
place  depends  to  a  large  extent  upon  the  activity  level  of  the  indivi- 
dual.  Reference  (5)  gives  typical  values  for  heat  production  of  man  for 
various  activities  as  shown  in  Table  I. 

Table  I 

Average  Metabolic  Rates 
for  Various  Activities 

Activity  Btu/hr,  Room  Temperature 

Between  60  F  and  90F 

Sleeping  280 

Lying,  relaxed  290 

At  rest,  seated  400 

At  rest,  standing  440 

Writing,  seated  430 

Typing,  seated  550 

Slow  movement  about  room  600 

Slow  walking  900 

Swimming  breat-stroke  1  mph  1620 

Wrestling  3100 
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Many  investigators  have  attempted  to  measure  effectively  the  rate 
of  production  of  body  heat.  DuBois  (6)  details  the  many  methods  used 
in  1939  the  principles  of  which  are  still  used  today. 

It  was  found  that  the  metabolic  rate  can  best  be  determined  by 

measuring  either  the  0  consumed  or  the  CO  produced.   The  first  method 

2  2 

is  the  one  generally  employed. 


Once  the  0  consumption  is  measured,  it  becomes  necessary  to 
evaluate  the  consumption  in  terms  of  units  of  energy  produced. 

The  caloric  value  of  a  liter  of  oxygen  consumed  by  an  individual 
will  depend  upon  his  previous  diet.   DuBois  (6)  points  out  that  a 

person  who  has  fasted  a  minimum  of  12  hours  will  have  a  respiratory 

CO2  produced  (Volume) 
quotient of  0.82. 

0  consumed  (Volume) 

Table  II 
Caloric  Value  of  1  Liter  0„ 


Caloric  Value 
Substance  1  Liter  0„ ,  K( 


Carbohydrate 

5.047 

Fat 

4.686 

Protein 

4.485 

For  Average  Fasting  R.Q. 

4.8 

Therefore,  employing  the  average  figure  of  4.8  for  the  caloric  value 
of  a  liter  of  0  can  give  results  whose  error  may  be  as  much  as  2% 
but  generally  less  than  1%  (6) . 


Radiative  Heat  Transfer 

The  net  radiative  heat  exchange  between  a  human  body  and  its 

surroundings  is  given  by  the  following  equation: 

R  =  o  A  1 (T  4  -  T   4)  |2 

r  1   +  A   (1  -  1)   s     mrt 
e,    A^   e 


1    Ae   e2 


where 


a  =  Stefan-Boltzman  constant  =  0.1713  x  10   Btu/ft2  hr  -  R4 

2 

A      =  Radiative  surface  area  of  the  body,  ft 

2 
A     =  Surface  area  of  the  environment,  ft 
e  ' 

e      =  Emissivity  of  the  human  skin 
1  ' 

e      =  Emissivity  of  the  surroundings 

Tg     =  Average  skin  temperature,  R 

Tmrt   =  Mean  radiant  temperature,  R 

However  when  the  surface  area  of  the  environment  is  large  compared 

to  the  area  of  the  person  the  equation  reduces  to: 

R  =  a  A   e   (T  4  -  T   *)  |3| 

r   1   s     mrt  '    ' 

The  radiative  surface  area,  A  ,  is  that  portion  of  the  total  sur- 
face area  which  is  seen  by  the  surroundings.   When  a  person  is  in  a 
sitting  position,  some  areas  of  his  skin  will  be  blocked  from  "seeing" 
the  surroudnings  by  other  areas  of  the  skin  and  only  a  part  of  the  total 
body  surface  area  will  be  contributing  to  the  radiative  heat  exchange. 

An  equation  for  determing  the  total  body  surface  area  of  an  indi- 
vidual based  upon  his  height  and  weight  was  developed  by  DuBois  and  DuBois  (7), 

A^     =   (0.108)  (W)0-425  (H)°-725 
b 

2 
A,      =  Total  body  surface  area,  ft 

W      =  Weight  of  an  individual,  lb 

H      =  Height  of  an  individual,  in 
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For  the  normal  seated  position,  Nielsen  and  Pedersen  (8)  found 
the  effective  radiation  area  to  be  approximately  60%  of  the  DuBois 
surface  area. 

Hardy  and  Muschenchin  (9)  have  shown  that  the  emissivity  of  human 
skin  approaches  that  of  a  black  body.   Since  the  surroundings  usually 
radiate  in  the  infrared  region  a  value  of  e  =  0.94  for  the  emissivity 
of  human  skin  can  be  used. 

Conductive  Heat  Transfer,  K 

Thermal  conduction  is  considered  as  a  transportation  of  heat 
energy  due  to  microscopic  motion,  such  as  molecular  motion  in  gases 
or  solids,  whereas  heat  transmission  connected  with  macroscopic  motion 
such  as  the  ordinary  flow  of  a  fluid  belongs  to  thermal  convection  (10). 

For  the  purposes  of  this  paper,  conductive  heat  transfer  was 
considered  to  be  that  associated  with  the  exchange  of  heat  between 
two  intimately  touching  solid  objects. 

Conductive  heat  loss  from  a  human's  body  would  then  occur  when  a 
portion  of  the  body  surface  area  was  in  contact  with  a  cooled  or  heated 
solid  surface.  A  person  seated  in  a  form  fitting  couch  would  have  a 
relatively  large  conductive  heat  loss.  The  minimum  conductive  loss 
would  probably  occur  during  swimming  or  sky-diving  when  the  body  touches 
no  solid  surfaces. 

When  a  person  is  seated  in  a  chair,  especially  one  with  webbing 
only  a  small  area  in  proportion  to  the  total  body  surface  area  is  in 
contact  with  the  chair.   For  such  a  small  area,  Jenssen  (11)  states 
that  the  conductive  coefficient  is  probably  not  too  different  from 
the  combined  convective  and  radiative  coefficients. 


10 
Therefore,  the  conductive  heat  loss  can  be  considered  negligible 
and  the  term  set  equal  to  zero  in  the  comfort  equation. 

Heat  Storage,  S 

The  storage  term,  "S",  may  be  plus  or  minus  as  the  temperature  of 

the  tissues  of  the  body  rise  or  fall.   To  find  the  amount  of  heat  gained 

or  lost  the  equation  S  =  m  C  dT  would  be  employed,  m,  the  mass  of  the 

P  dt 
human  body,  is  determined  by  weighing  the  subject.   C  ,  the  specific 

heat  of  the  human  body  tissues  is  usually  taken  to  be  0.83  Btu/lbs  F  (12)  (13) 

The  average  temperature  of  the  human  body  tissues  has  been  given 

as 

T    =  0.65  T       +  0.35  T     by  Burton  (14). 
avg        rectal        skin 

Another  formula,  only  slightly  different,  was  given  as 

T    =  0.8  T  +  0.2  T  by  DuBois  (6). 
avg        r        s  J 

Winslow  and  Herrington  (15)  compared  the  two  equations  and  state 
that  for  the  ambient  temperature  range  77-93.5  F,  the  DuBois  formula 
is  better  but  for  the  lower  ambient  temperature,  the  Burton  formula 
is  better. 

With  a  little  reflection,  it  can  be  seen  that  the  storage  term 
could  be  positive  or  negative  for  a  short  period  of  time  and  the  human 
body  remain  in  the  comfort  zone.  However,  for  longer  periods  of  time, 
if  the  body  continues  to  lose  or  gain  heat,  it  will  leave  the  zone  of 
thermal  equilibrium  and  the  comfort  equation  will  no  longer  be  applicable. 

Since  equation  |l  |  is  to  be  used  in  cases  where  thermal  equili- 
brium is  to  be  maintained  the  storage  term  S  will  be  set  equal  to  zero. 

Evaporative  Heat  Loss,  E 

The  loss  of  heat  by  evaporation  of  moisture  from  the  surfaces  of 
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the  skin  and  lungs  is  one  of  the  easiest  terms  in  the  heat  balance  to 
determine,  yet  the  most  difficult  to  describe  analytically. 

By  examining  experimental  data  Gagge  (16)  introduced  the  following 
equation  for  calculating  the  heat  loss  due  to  evaporation. 

E  =   (Wy)  A.  (P  -  RH  x  P  )CF  |4  | 

where 

W  =  Fraction  of  body  area  completely  wet 

y  =   Coefficient  of  heat  transfer  by  evaporation,  containing  the 
constants  for  evaporization,  air  motion  and  direction 

Btu/ft  -hr-psia  vapor  pressure 

o 
A^  =  Total  body  area  (DuBois  formula) ,  ft 

P  =  Saturation  vapor  pressure  at  skin  temperature,  psia 

RH  =  Relative  humidity,  % 

Pa  =  Saturation  vapor  pressure  at  ambient  temperature,  psia 

CF  =  Conversion  factor,  (1.90633)   Btu  cm  Hg  m2 

KCal  psia  ft2 

Gagge  points  out  that  the  (Wy)  term  cannot  be  separated  into  its 
components.  Winslow  and  Herrington  (15)  clarify  this  by  saying  that 
this  index  concerns  the  physiological  reaction  of  the  body  governing 
the  rate  of  sweat  secretion  and  not  merely  what  area  of  the  body  is 
covered  with  sweat. 

Gagge  found  there  is  a  maximum  value  for  (Wy)  of  29.9  K  ,  ,  „     2 
for  the  total  moist  surfaces  of  the  body,  beyond  which  the  evaporative 
heat  loss  does  not  increase.   As  can  be  seen  by  examining  the  above 
equation,  the  maximum  evaporative  heat  loss,  when  (Wy)  is  maximum,  will 
depend  upon  the  skin  temperature,  air  temperature,  and  relative  humidity 
for  low  air  velocity. 
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Gagge  also  found  that  (Wp)  does  not  go  to  zero,  but  remains  at 

approximately  12-15%  of  (Wy)  maximum  regardless  of  how  cold  the  skin 

2 
becomes.   A  value  of  (Wy)  =  4.7  KCal/cmHg-m  can  be  used  to  represent 

the  minimum  value.  The  resulting  minimum  evaporative  heat  loss  is 

dua  to  the  evaporation  of  water  which  diffuses  through  the  epidermis 

from  the  deeper  layers  of  skin  and  the  evaporation  of  moisture  from 

the  moist  surfaces  of  the  respiratory  system.   Since  vapor  pressure  is 

the  driving  force  for  this  minimum  evaporative  cooling,  it  will  be 

independent  of  the  total  ambient  pressure,  and  of  the  gas  surrounding 

the  skin. 

Convective  Heat  Transfer,  C 

Heat  transfer  by  convection  between  a  solid  boundary  (skin)  and 
a  fluid,  liquid  or  gas,  is  evaluated  by  means  of  the  equation: 

C  -  Achc  <Ts  -  Ta>  '5| 

where 

C  =  rate  of  convective  heat  loss,  Btu/hr 

Ac  =  convective  area,  ft2  =  Av 

Ts  =  average  skin  temperature,  F 

Ta  ■  dry  bulb  temperature  of  ambient  fluid,  F 

hc  =  average  convective  heat  transfer  coefficient,  Btu/hr  ft2  F 
The  convection  equation  in  this  form  seems  quite  simple.   The  simplicity 
is  misleading,  however,  due  to  the  difficulty  in  evaluating  hc.   The 
convective  heat  transfer  coefficient  is  actually  a  complicated  function 
of  the  thermal  properties  of  the  fluid,  the  geometry  of  the  system  and 
the  fluid  flow.  This  coefficient  may  be  expressed  as  a  function  of 
the  following  parameters: 
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h    =  f(L,  v,   p,  k,  C  ,  AT,  3,  g,  V) 
c  p 

where 

L  =   characteristic  geometrical  dimension,  ft 

|i 

p  =  density  of  the  fluid,  lb/ft" 

k  =  thermal  conductivity  of  the  fluid,  Btu/hr-ft-F 


=  viscosity  of  the  fluid,  lb/hr-ft 

3 


C   =  specific  heat  at  constant  pressure  of  the  fluid,  Btu/lb-F 
P 

AT  =  temperature  difference  between  body  surface  and  ambient  fluid,  F 

8   =  coefficient  of  volumetric  expansion,  F 

V   =  velocity  of  fluid,  ft/hr 

g   =  acceleration  due  to  gravity,  ft/hr2 

In  order  to  facilitate  handling  of  so  many  properties  of  systems, 
scientists  through  the  years  have  developed  sets  of  dimensionless 
relationships.   One  such  combination  called  the  Nusselt  modulus  or 
Nusselt  number,  Nu,  combines  the  convective  heat  transfer  coefficient 
hc,  the  significant  dimension  L,  and  the  thermal  conductivity  of  the 
fluid,  K. 

Nu  =  hc  L  1 6 1 

Other  dimensionless  groups  with  wich  the  paper  will  be  concerned 


are: 


Prandtl  number  (Pr)   =  C  y 

P 


Grashof  number  (Gr)   =  p2gL^  (AT) 

U2 


Reynolds  number  (Re)   =  L_Vp 

V 
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The  complicated  geometrical  shape  of  the  human  body  makes  it 
extremely  difficult,  if  not  impossible,  to  determine  the  value  of  the 
convective  heat  transfer  coefficient  accurately.   Some  attempts  have 
been  made  (17) ,  however  the  results  are  valid  only  in  special  cases 
and  even  then  are  not  conclusive. 

The  next  step  is  to  approximate  the  human  body  by  a  simple  geome- 
trical shape,  for  which  correlation  groupings  of  the  appropriate 
dimensionless  groups  and  the  resulting  equations  exist.   The  two  shapes 
considered  were  spheres  and  cylinders  with  cylinders  generally  being 
used. 

Before  the  actual  equation  is  written,  the  type  of  fluid  flow 
should  be  considered. 

Convection  heat  transfer  can  be  divided  into  two  types;  forced 
convection  and  free  or  natural  convection,  depending  upon  the  velocity 
of  the  fluid. 

Forced  Convection 

Forced  convection  is  usually  considered  to  be  that  which  occurs 
when  the  velocity  of  the  fluid  is  great  enough  to  have  an  effect  upon 
the  heat  transfer  process.   A  person  outdoors  on  a  windy  day  would 
lose  heat  by  the  process  of  forced  convection.   McNall  and  Sutton  (18) 
suggest  that  for  air  velocities  of  50  ft/min  or  greater,  an  equation 
for  forced  convection  should  be  used. 

Assuming  that  the  fluid  flows  in  a  horizontal  direction  normal 
to  the  cylinder,  McAdams  (19)  suggests  an  equation  of  the  type 

Nu  =   4  (Re)n  (Pr)m  |7| 
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where 

<j>   =  denotes  a  functional  relationship  and  n  and  m  are  constants. 

Free  Convection 

Free  convection  heat  transfer  occurs  whenever  a  body  is  placed 
in  a  fluid  at  a  higher  or  lower  temperature  than  that  of  the  body. 
As  a  result  of  the  temperature  difference  heat  flows  between  the  fluid 
and  the  body  causing  a  change  in  the  density  of  the  fluid  layers  in 
the  vicinity  of  the  surface.   The  difference  in  density  leads  to  a 
downward  flow  of  the  heavier  fluid  and  an  upward  flow  of  the  lighter 
fluid.   If  the  motion  of  the  fluid  is  caused  solely  by  differences  in 
density  resulting  from  temperature  gradients,  without  the  aid  of  a 
pump  or  fan,  the  associated  heat  transfer  mechanism  is  called  natural 
or  free  convection.   Free  convection  currents  transfer  internal  energy 
stored  in  the  fluid  in  essentially  the  same  manner  as  forced  convection 
currents.  However,  the  intensity  of  the  mixing  motion  is  generally 
less  in  free  convection,  and  consequently  the  heat  transfer  coefficients 
are  lower  than  in  forced  convection  (20) . 

In  many  practical  situations,  man  is  seated  and  relatively  motion- 
less while  at  work.   In  such  cases,  free  convection  equations  usually 
will  be  applicable. 

For  free  convection  heat  transfer  McAdams  (19)  suggests  an  equation 

of  the  form: 

0  25 
Nu      =     <j>(Gr  Pr)    *      ,  which  for  a  vertical  cylinder  becomes 

Nu      =     0.59    (Gr  Pr)0*25  1 8 1 
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When  the  buoyancy  is  the  only  driving  force  the  fluid  velocity 
is  determined  entirely  by  the  quantities  contained  in  the  Grashof 
number.   Therefore,  the  Reynolds  number  is  superfluous  for  free  con- 
vection and  is  dropped  from  the  equation. 

For  a  horizontal  cylinder  McAdams  (19)  recommends  Nu  ■  0.53 

0  25 
(Gr  Pr)   '  .  When  an  individual  is  lying  flat  or  standing  upright, 

the  choice  is  obvious.   However,  for  the  seated  person  the  choice  is 

more  complicated  since  the  thighs,  feet,  hands  and  arms  are  nearly 

horizontal  while  the  rest  of  the  body  can  be  considered  vertical.   Since 

the  horizontal  surface  area  would  be  approximately  45%  of  the  total 

surface  area  the  equation,  Nu  =  { (0.45) (0.53)  +  (0.55) (0.59) }  (Gr  Pr)  * 

resulting  in  Nu  =  0.56  (Gr  Pr)  *  ->.   This  is  very  close  to  the  value 

of  Nu  =  0.555  (Gr  Pr)0,25  suggested  by  Kreith  (20)  for  vertical  cy- 

g 

linders  when  the  Grashof  number  is  less  than  10  as  it  is  in  experiments 

involving  the  human  body  in  thermal  equilibrium. 

The  convection  heat  transfer  equation  now  reads: 


c 

Ab  hc  AT 

where 

hc    = 

k  (40  £ 

L 

by  defining 

Cc    = 

hc 
AT0.25 

V 


|9 


.3  8  B  P2  AT  Pr)0'25 

V1  llOJ 


1  25 

Therefore,  C  =  A^  C   (AT)  '  |ll| 

The  (AT)  term  has  been  removed  from  the  convection  coefficient  and 
a  more  general  form  has  been  achieved,  namely 

C     -  K   oc  <♦>   (g  3  P2  Pr)0'25 

cc    -  10.25      5 I   I 

Vz  .  .  12 
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The  terms  in  C  which  are  temperature  dependent  are  evaluated  at 
film  temperature  T  which  is  the  average  of  the  temperature  of  the 
skin  and  the  temperature  of  the  ambient  fluid.  Those  terms  in  C 
which  are  independent  of  the  fluid  considered  are: 

L    =  unknown,  ft 

(j>    =  unknown,  dimensionless 

g    =  4.17  x  10  ,  ft/hr2 

3    =  1/T  ,  1/R 

For  air  all  properties  were  obtained  from  reference  (21). 

For  helium 

Pr   =  0.6871  and  is  a  constant  for  the  temperatures  and 
pressures  considered  (22) . 

In  the  following  relationships  for  p,  u,  and  k  as  function  of 

temperature  and  pressure,  T  =  temperature  in  degrees  R,  and  P  ■ 

pressure  in  psia. 


2.6829  Tf  +  P  3°'25  (0.3059  -  1.824  B°*5  -  0.8223) 


lb/ft3  (23) 


y    =  8.315  (Tf)0,647  10  A  lb/hr  ft  (24) 

=  0.0792  +  0.962  10~4  (T,  -  459.59)  +  10  10"5   (  p    -  1)  (23) 

14.696 

This  leaves  only  <J>  and  L  to  be  defined.   Since  they  cannot  be 

specifically  determined,  the  combination   <f>    needs  to  be  found. 

Nielsen  and  Pedersen  (8)  found  by  experimentation  a  value  for 
C  for  the  human  body  in  a  seated  position.   The  experiments  were  run 
with  an  air  velocity  less  than  20  ft/min;  therefore,  the  equation  for 
free  convection  heat  transfer  would  apply.   The  value  they  found  for  C 


c 


in  metric  units  is  2.03  cal      .  When  this  is  converted  to  English 

m2hr  C1,25 
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units  it  b( 

=comes  0.359  Btu     ,  25* 
ft7  hr  F  * 

By  se 

tting  this  equal  to  equation  |l2'|  for  C  with  the  values  for 

c 

the  thermal  properties  of  helium  at  90F  and  one  atmosphere    C  He  -  0.754. 

Final  Heat 

Balance  Equation 

Substituting  equations  |4|  | 11 |  and  |3|  for  their  respective 

symbols  into  equation  |l|  ,  the  following  relationship  is  obtained. 

M  = 

(Wy)A   (P  -  RH  P  )  CF  ±  A  C  (T  -  T  )1*25  ±  aA  e,  (T  A-T   4)   1 13 1 
Pbs       a       bCs    a         rls   mrt    '   ' 

where 

(Wy) 

=  Constant  representing  the  fraction  of  body  area  completely 

wet  and  containing  the  constant  for  evaporation.  (Wy)  minimum  = 
4.7  KCal/cmHgm2 

\ 

2 
=  Total  body  surface  area,  ft 

P 
s 

=  Saturation  vapor  pressure  at  skin  temperature,  psia 

RH 

=  Relative  humidity,  % 

P 
a 

=  Saturation  vapor  pressure  at  ambient  temperature,  psia 

CF 

2 
=  Conversion  factor,  (1.90633)  Btu  cmHgm 

KCal  psia  ft 

2    1  25 
Cc     =  Modified  coefficient  of  convective  heat  transfer  Btu/ft^hr  F  * 

T      =  Skin  temperature 

T      =  Dry  bulb  temperature  of  the  ambient  gas,  F 

3. 

-8       2    4 
a  =  Stefan-Boltzman  constant  0.1713  x  10   Btu/ft  hr  R 

2 

A      =  Radiative  surface  area  of  the  body,  ft 

c,  =     Emissivity  of  the  human  skin,  0.94 

T      =  Mean  radiant  temperature,  R 

One  purpose  for  deriving  equation  |13|  is  to  be  able  to  predict  the 
environmental  conditions  for  thermal  equilibrium  and  hopefully  comfort. 
It  should  be  noted  that  equation  |l3|  applies  only  when  the  person  is 
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sedentary  or  near  sedentary  and  free  convection  terms  are  applicable. 
In  deep  diving  expeditions,  it  is  assumed  the  occupants  will  be  relatively 
inactive  while  in  the  diving  capsule  and  will  save  their  energy  for 
their  work  in  the  water. 

In  order  to  apply  equation  | 13 |  an  average  value  for  skin  tempera- 
ture must  be  assumed.   Although  man  may  be  in  thermal  equilibrium  over 
a  range  of  skin  temperatures,  91-93  F  was  found  to  be  optimum  (24). 
Once  a  skin  temperature  has  been  chosen,  the  heat  balance  equation 
can  be  solved  to  determine  the  environmental  dry  bulb  temperature  necessary 
to  produce  thermal  comfort  for  a  given  mean  radiant  temperature. 


20 

EXPERIMENTAL  METHOD 

Six  tests  were  conducted  involving  a  seated  subject  exposed  to  a 
helium  environment  contained  in  a  polyethylene  test  cell  using  the 
equipment  and  techniques  described  in  the  body  of  this  report.   In 
addition,  a  test  was  conducted  in  the  test  cell  with  air  as  the  ambient 
gas  and  two  tests  were  conducted  in  air  using  the  same  equipment  but 
outside  the  polyethylene  test  cell.   The  subject  recorded  his  sub- 
jective thermal  comfort  during  the  tests.   Also  the  parameters  necessary 
to  determine  the  subject's  heat  balance  were  measured  and  calculations 
according  to  theory  were  made. 

Subiect 

The  subject  was  an  engineering  graduate  student  25  years  of  age  in 
good  health  and  good  physical  condition. 

His  build  was  average  with  a  height  of  70.25  inches  and  a  weight 
of  157  pounds.   While  the  subject  would  not  be  classified  as  a  trained 
subject,  he  was  not  naive  to  environmental  and  physiological  experi- 
ments and  their  procedure. 

In  order  to  use  the  caloric  value  of  oxygen  determined  for  a  fasting 
person  (see  p.  7),  the  subject  did  not  eat  after  8:00  p.m.  of  the 
evening  preceding  the  test.   The  subject  had  a  minimum  of  six  hours 
sleep  the  preceding  night  and  refrained  from  strenuous  activity  the 
evening  before  and  the  morning  of  the  test.   In  the  weeks  preceding 
the  tests,  the  subject  followed  the  pre-test  fasting  procedure  on 
selected  days  and  his  metabolic  rate  was  determined  for  the  sedentary 
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sitting  position  in  air  at  a  mean  temperature  of  74  F.  This  acquainted 
the  subject  with  the  experimental  procedure. 

Equipment 

The  experiment  was  conducted  in  the  Environmental  Test  Room  of 
the  Environmental  Research  Institute,  Kansas  State  University.   The 
room  is  12  x  24  x  8  feet  in  its  present  configuration.  The  walls, 
floor  and  ceiling  consist  of  aluminum  panels  backed  with  copper  tubing. 
By  controlling  the  temperature  of  the  liquid  flowing  through  the 
tubing,  the  surface  temperature  of  the  enclosure  can  be  maintained 
within  ±1  F  dry  bulb  and  wet  bulb. 

A  polyethylene  bag  (Figure  1)   7  x  5.5  x  4.5  feet  was  constructed 
by  heat  sealing  "Visqueen",  a  product  of  Ethyl  Corporation,  Baton  Rouge,  La. 

The  bag  will  hereinafter  be  referred  to  as  the  test  cell. 

At  first  the  test  cell  was  without  openings  except  for  a  slit 
approximately  4  feet  wide,  1  foot  above  the  floor.  This  opening 
allowed  the  subject  to  enter  the  test  cell. 

Since  the  dew-point  of  the  gas  in  the  test  cell  was  to  be  con- 
trolled by  maintaining  the  floor  of  the  test  cell  at  a  low  temperature, 
it  was  necessary  that  the  floor  temperature  be  uniform  and  that  good 
heat  transfer  exist  between  the  test  room  floor  and  test  cell  floor. 
With  the  polyethylene  bottom  taped  to  the  test  room  floor,  this  could 
not  be  done.   A  4  x  8  ft  sheet  of  1/8  inch  Alcoa  Aluminum  was  then 
placed  on  the  test  room  floor,  the  bottom  cut  out  of  the  test  cell, 
and  the  edges  of  the  test  cell  taped  to  the  aluminum  with  polyethylene 
tape  and  ducting  tape.  The  aluminum  sheet  temperature  was  within 
±0.5  F  of  the  test  room  floor  temperature  and  could  be  controlled  easily. 
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Figure  1.      Inflated  Polyethylene  Test  Cell  Supporting  Structure 
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The  supporting  structure  shown  in  Figure  1  was  used  to  support 
the  test  cell  between  tests  and  during  the  inflation  process.  The 
test  cell  remained  attached  to  the  structure  during  the  actual  test 
although  it  was  probably  not  necessary  due  to  the  buoyancy  of  the 
helium.   A  Fairbanks-Morse  beam  balance  scale,  accurate  to  ±0.01  lbs, 
was  placed  inside  the  test  cell  and  a  lawn  chair  with  most  of  the 
webbing  removed  placed  upon  it  and  taped  in  place  (see  Figure  2). 

Twenty-three  thermisters,  Model  409,  manufactured  by  the  Yellow 
Spring  Instrument  Company,  Yellow  Springs,  Ohio,  were  used  to  monitor 
skin  temperature  and  the  temperature  of  the  test  cell.   And  one  YSI 
Model  410  thermister  was  used  to  monitor  deep  body  temperature.   Since 
the  thermisters  were  equipped  with  only  six  foot  leads,  a  junction 
board  was  built  and  attached  to  the  scale  (see  Figure  2).   The 
junction  board  leads  were  passed  through  the  umbilical  connector  to 
the  outside  where  they  were  plugged  into  the  read-out  device. 

A  United  Systems  Corporation,  Dayton,  Ohio,  Series  500  Digitec 
Thermometer  was  used  to  read  the  thermister  temperatures  and  a  Digitec 
Series  500  Digital  Printer  was  used  to  simultaneously  record  the 
readings  (see  Figure  3).   The  accuracy  of  the  thermometer  was  ±0.15  F. 

The  calibration  of  the  thermisters  was  checked  by  placing  them 
with  a  mercury  thermometer  in  a  water  vessel  and  allowing  them  to 
reach  equilibrium.   The  thermister  readings  checked  with  each  other 
and  the  mercury  thermometer  to  within  ±0.1  F. 

A  face  mask  Model  133-054  adult  size  obtained  from  the  Puritan 
Compressed  Gas  Corporation,  Kansas  City,  Missouri  was  worn  by  the 
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Figure  2.   Subject  Seated  in  Test  Configuration  Showing  the  Mask,  Valve 
Humidity  Sensor,  Scale,  and  The mister  Junction  Board,  Test 
Cell  Removed. 
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subject  during  all  test  runs  (see  Figure  2).  The  inflatable  rim 
allows  the  straps  to  pull  the  mask  tight  enough  to  seal,  yet  not 
be  uncomfortable. 

The  mask  was  connected  with  practically  no  dead  space  to  a 
modified  Otis-McKerrow  two-way  intake-exhaust  valve  (see  Figure  2). 
Each  side  of  the  valve  was  connected  to  the  air  outside  the  bag  by  a 
1%  inch  flexible  vinyl  tube  which  passed  through  the  umbilical  connec- 
tion. 

A  chain-compensated  172  liter  gasometer  manufactured  by  the 
Collins  Corporation  was  used  to  measure  oxygen  consumption.   The  level 
of  oxygen  in  the  gasometer  tank  was  transformed  to  an  output  voltage 
which  was  recorded  on  an  Omnigraph  Corporation,  Bellaire,  Texas, 
Model  HR-95  X-Y  Recorder. 

Medical  oxygen  was  used  in  the  gasometer  while  medical  helium 
provided  the  environment  for  the  experiment.  The  equipment  needed  by 
the  subject  consisted  of  the  following  items:   pencil,  paper,  foam 
cushions,  extra  scale  weights  and  relative  humidity  sensor.   The  pencil 
and  paper  were  needed  for  taking  notes  and  recording  data,  the  foam 
cushions  were  placed  over  sharp  portions  of  the  scale  while  the  bag 
was  collapsed  and  the  extra  weights  were  used  to  keep  the  subject's 
weight  at  a  value  that  could  be  easily  read. 

Pre-Test  Preparation 

The  subject  entered  the  pre-test  room  which  was  maintained  at  a 
temperature  between  72-76  F,  at  least  thirty  minutes  prior  to  the  time  . 
for  the  test.  At  this  time  he  changed  clothes  putting  on  the  swimsuit 
he  would  wear  during  the  test. 
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For  the  first  three  tests  in  the  series,  the  thermister  harness 
was  taped  on  while  the  subject  was  in  the  pre-test  room.   Later,  after 
a  few  of  the  thermisters  had  been  moved  from  the  body  to  locations  on 
the  chair  and  exposed  to  the  gas  in  order  to  determine  the  thermal 
gradients  within  the  test  cell,  the  thermisters  were  taped  on  after 
the  subject  was  seated  in  the  chair  within  the  test  cell.   It  was  found 
that  this  required  about  the  same  amount  of  time  and  eliminated  the 
re-taping  of  those  thermisters  which  were  pulled  loose  while  the 
subject  entered  the  test  cell  and  was  seated  in  the  chair.  The  rectal 
probe  had  been  inserted  by  the  subject  while  dressing. 

After  the  subject  was  seated  and  all  thermisters  checked,  the  air 
bubble  on  the  scale  was  checked  to  insure  that  the  scale  was  level, 
and  if  not,  the  scale  was  leveled.   The  next  step  was  to  put  on  the 
breathing  mask  and  to  check  it  for  leaks.  This  was  done  by  the  subject 
by  placing  a  hand  on  the  exhaust  port  of  the  valve  and  blowing.  The 
hoses  were  then  connected  to  the  mask  and  the  assistant  left  the  test  cell. 

At  this  point  the  test  cell  was  unhooked  from  its  support  and 
allowed  to  collapse  around  the  subject  and  scale.   That  section  of  the 
test  cell  containing  the  opening  was  flattened  against  the  floor  and 
the  opening  sealed  with  polyethylene  tape. 

After  the  cushions  had  been  properly  placed  on  the  scale  to  avoid 
puncturing  the  test  cell,  a  vacuum  cleaner  was  attached  to  a  tube  in- 
cluded in  the  umbilical  for  this  purpose,  and  the  suction  used  to  draw 
the  air  from  the  test  cell.  As  shown  in  Figure  4,  the  test  cell, 
after  the  vacuum  cleaner  was  turned  off,  was  very  tight  around  the 
subject.   Helium  was  then  leaked  into  the  test  cell  at  a  slow  rate 
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Figure  4.   Polyethylene  Test  Cell  Collapsed  Showing  the  Degree  of 
Evacuation  of  Air  From  the  Test  Cell 
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until  the  test  cell  had  moved  several  inches  from  the  subject.   The 
vacuum  cleaner  was  turned  on  again  and  the  test  cell  evacuated  for 
a  second  time.  The  vacuum  cleaner  was  turned  off  and  the  helium  flow 
started  once  again.   The  helium  remained  flowing  continuously  at  various 
rates  until  the  end  of  the  test. 

This  evacuation  process  and  the  fact  that  a  constant  positive 
pressure  was  maintained  inside  the  test  cell  were  the  bases  for  the 
assumption  that  the  environment  was  for  practical  purposes  100%  helium. 

The  test  cell  was  then  slowly  filled  with  helium;  this  required 
approximately  15-20  minutes. 

Testing  Procedure 

The  actual  test  procedure  began  as  the  test  cell  was  filling  slowly 
to  its  inflated  configuration,  (see  Figure  1).   During  this  time,  the 
subject's  weight  was  taken  and  an  initial  relative  humidity  and  tem- 
perature check  made  by  the  subject.   The  assistant  monitoring  the  skin 
thermister  read-out  and  printer  noted  the  various  skin  temperatures  and 
compared  them  with  the  previous  tests.   If  any  were  found  to  be  approximately 
2  degrees  F  high  or  low,  a  check  was  made  to  see  if  the  sensors  were 
still  in  contact  with  the  skin,  between  the  skin  and  chair  or  not  in  the 
proper  location.   Then,  the  skin  and  gas  temperatures  were  allowed  to 
reach  equilibrium. 

The  test  room  walls,  ceiling,  and  floor  were  set  at  a  specified 
temperature  before  testing  began  and  remained  constant  during  the  test. 
The  temperature  of  the  air  in  the  test  room  was  varied  as  needed  to  help 
control  the  test  cell  gas  temperature.   The  only  other  control  used  was 
the  rate  of  helium  flow  into  the  test  cell.   The  expansion  of  the  helium 
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proved  to  be  a  very  good  heat  sink  and  was  used  to  control  the  gas 
temperature. 

Once  the  gas  and  skin  temperature  were  stable,  the  actual  test  began. 
The  subject's  initial  weight  was  determined  and  the  temperature  and  rela- 
tive; humidity  racordad.   Temperatures  were  recorded  every  ten  minutes; 
subject  weight,  relative  humidity  and  temperature  of  the  gas  were  recorded 
every  twenty  minutes. 

The  subject  quickly  learned  that  he  must  assume  the  same  position 
each  time  he  recorded  his  weight  since  the  tubes  connected  to  the  valve 
were  taped  to  the  scale.   The  position  of  the  head  was  very  important 
and  was  maintained  constant  by  resting  the  chin  on  the  left  hand  and 
sighting  across  the  face  mask  and  valve  to  the  same  point  each  time. 
Also,  it  was  learned  that  the  scale  was  so  sensitive  that  all  persona 
in  the  test  room  had  to  refrain  from  walking  during  weighing. 

Once  each  test,  a  metabolic  rate  was  determined  by  measuring  the 
subject's  oxygen  consumption.   Suffice  it  to  say  that  the  only  differ- 
ence to  the  subject  when  the  oxygen  consumption  was  being  measured  was 
the  difference  in  humidity  between  the  saturated  oxygen  and  the  room 
air  usually  inhaled  by  the  subject  and  that  the  subject  refrained  from 
taking  any  measurements  in  order  to  obtain  a  more  regular  oxygen  con- 
sumption curve. 

As  the  wearing  of  the  mask  was  not  a  new  experience  to  the  sub- 
ject, its  affect  on  the  metabolic  rate  should  have  been  constant 
throughout  the  test  series. 

After  the  metabolic  rate  was  determined,  the  test  proceded  until  a 
weight  loss  of  at  least  0.1  pounds  was  recorded.   This  was  usually  a 
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period  of  two  hours  after  stabilization. 

At  the  conclusion  of  the  test,  the  helium  was  turned  off,  the  tape 
was  removed  from  the  test  cell  opening,  and  the  vacuum  cleaner  used  to 
draw  air  through  the  test  cell.   In  this  way,  the  helium  environment 
was  diluted  and  the  subject  could  remove  the  mask  before  leaving  the 
chair. 

Evaluation  of  Physical  Parameters 

In  order  to  determine  the  radiation  term  of  equation  | 13 |  ,  the 
following  parameters  were  to  be  evaluated:   skin  temperature,  the  mean 
radiant  temperature  of  the  surroundings  and  the  radiation  area  of  the 
subject. 

The  subject  in  the  experiment  weighed  147  pounds  with  a  height 
of  70.25  inches.   Using  the  DuBois  formula  his  total  body  surface  area 
was  calculated  to  be  20.3  square  feet.   Therefore  A  =  (0.60) (20.32)  = 
12.2  square  feet. 

Since  all  parts  of  the  human  body  are  not  the  same  temperature,  it 
is  necessary  to  calculate  an  average  skin  temperature.   This  was  done 
by  measuring  the  skin  temperature  at  nineteen  points  after  Hardy  and 
DuBois  (25),  (see  Figures  5  and  6).  The  skin  temperatures  for  the 
various  parts  are  then  weighted  according  to  the  following  table  to 
give  the  mean  skin  temperature.   The  weighting  factor  corresponds  to 
the  percent  of  total  body  area  for  that  part. 

The  following  equation  was  used  to  compute  the  skin  temperature. 
Tg  =  0.07  T±  +   0.14  T2  +  0.05  T3  +  0.07  T4  +  0.13  T5  +  0.19  Tfi  +  0.35  I   | 14 
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Table  III 
Skin  Temperature  Weighting  Factors 
Part  of  the  body       Temperature     Weighting  Factor 

0.07 
0.14 
0.05 
0.07 
0.13 
0.19 
0.35 
Total    1.00 

An  equation  for  mean  radiant  temperature  is  given  by  Raber  and 
Hutchinson  (26)  as 

Tmrt  "  F12  T2  +  F13T3+  '"  Fn  Tn  1 15 1 

■ 

where 

F     =  the  shape  factor  from  body  one,  in  this  case  the  seated 
•k3,      person  ,  to  body  n 

T     =  Temperature  of  body  n,  F 

This  equation  for  the  evaluation  of  T   .  is  valid  when  the  surfacing 

n  mrt 

materials  have  an  emissivity  in  excess  of  0.90. 

Since  the  floor  was  maintained  sometimes  at  a  lower  temperature  than 
the  walls  and  ceiling,  it  was  necessary  to  evaluate  the  shape  factor  from 
the  seated  test  subject  to  the  floor.   This  was  done  using  the  method  of 
Tripp  (27).   The  man's  radiative  surface  area  was  represented  by  a 
cylinder  4  feet  high  with  a  radius  of  0.465  feet.   This  cylinder  placed 
in  the  location  shown  by  Figure  7  approximates  the  location  of  the  test 
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Figure  5.     Subject  in  Thermister  Harness,  Front  View 
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Figure  6.     Subject  in  Thermister  Harness,  Back  View 
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subject  (A  feet  from  the  west  wall  and  8  feet  from  the  south  wall). 

2*' 


Figure  7.   Test  Room  Floor  Showing  Location  of  Subject 
and  Division  of  Floor  for  Determination  of 
shape  factors. 
In  order  to  use  the  graphs  provided  by  Tripp,  the  test  room  floor 

was  divided  into  the  eight  sections  as  shown  in  Figure  7.   The  shape 

factor  was  then  found  from  the  right  circular  cylinder  to  each  of  the 

right  triangles.  The  summation  of  these  shape  factors  yielded  a  shape 

factor  from  the  man  to  the  floor  of  F    ....    =0.33.   The  shape  fac- 

man-floor 

tor  from  the  man  to  the  walls  and  ceiling  then  equals  0.67. 

The  walls  of  the  environmental  test  room  are  covered  with  at 
least  five  coats  of  light  green  enamel  (glossy)  the  emissivity  of  which 
is  approximately  equal  to  0.95.   The  floor  is  covered  with  a  vinyl 
carpet  for  which  the  emissivity  is  also  assumed  to  be  0.95  (28). 

A  spectral  analysis  of  the  polyethylene  used  in  the  construction 
of  the  test  cell  was  done  by  Mr.  H.  T.  Leventhal  (29).   This  analysis 
showed  the  test  cell  to  have  a  transmittance  of  0.84  in  the  infrared 
range.   The  film  was,  therefore,  considered  to  be  "transparent"  for 
purposes  of  radiation  heat  transfer. 

Ag  =  1152  square  feet  for  the  test  room  (12  x  24  x  8  feet) ,  there- 
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fore  the  term    1 becomes  1 =     1 -  e^ 

1  +  A  ,1  1     +   12.2  (0.05)   1  +0.0005 

e±         r  %  -V  tx       1152         ex 

(within  the  accuracy  of  this  report)  and 

R  =  a   A  e,  (T  A  -  T  M. 

r  Is     mrt 

Due  to  the  low  floor  temperature  and  low  air  circulation,  a  thermal 
gradient  existed  in  the  ambient  gas.   The  gas  temperature  was  measured 
with  skin  thermisters  attached  to  the  chair  and  exposed  to  the  atmosphere; 
also,  a  thermister  was  attached  to  the  rod  shown  extending  behind  the 
chair  in  Figure  2.  The  arithmetic  mean  of  the  measurements  was  used 
as  T  .  The  AT  over  the  approximate  height  of  the  subject  varied  from 
3.8  to  0.8  F. 

In  an  attempt  to  simulate  the  minimum  condition  of  heat  transfer 
by  conduction  the  chair  used  in  the  experiment  was  modified  such  that 
only  a  small  area  was  in  contact  with  any  part  of  it.  Therefore,  the 
conductive  heat  loss  K  was  set  equal  to  zero  in  the  heat  balance  equa- 
tion. 

The  experiment  was  conducted  such  that  the  mean  body  temperature 
did  not  change  appreciably  and  the  storage  term  S  was  set  equal  to  zero. 
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RESULTS  AND  DISCUSSION 

Theoretical  Analysis 

The  heat  balance  equation,  | 13 |  p.  18  can  be  used  to  describe  the 

experiment  of  this  report  by  substituting  values  of  A^  =  20.3  square 

feat  and  A  -  12.2  square  feet  which  were  calculated  for  the  subject. 

For  each  selected  value  of  Tgl  a  value  of  T^  was  calculated  for  all 

combinations  of  T  =  60,  65  ••'110  F,  RH  =  20,  40,  60,  80%  and  (Wu)  =  4.7 
a 

and  29.9  KCal/cm  Hg-m2. 

The  solutions  of  equation  1 13  j  for  air  and  helium  at  one  atmos- 
phere pressure  and  for  helium  at  20  atmospheres  pressure  with  assumed 
skin  temperatures  of  90,  92,  94  and  96  F  are  shown  in  Figures  9  through 
17.   These  data  were  generated  by  an  IBM  1410  digital  computer.   A 
copy  of  the  Fortran  program  is  included  in  the  appendix. 

Since  most  sedentary  people  are  comfortable  when  their  only  per- 
spiration is  of  the  insensible  type,  a  value  of  (Wm  )  Minimum  has  been 
chosen  to  represent  the  minimum  wetted  area.   This  indicates  a  minimum 
evaporation  rate  which  is  shown  by  the  shaded  band  of  Minimum  in 
Figures  9  through  17  and  is  therefore  called  the  "comfort  zone"  for 
the  designated  skin  temperature.   The  region  below  this  band  is  the 
zone  of  body  cooling,  and  although  a  person  may  "live"  at  this  point, 
he  is  compensating  for  his  heat  loss  by  means  other  than  normal  vaso- 
motor regulation  against  cold.  The  area  above  the  comfort  band  is 
the  area  of  evaporative  regulation.   In  this  area,  the  sweat  rate  is 
increased  and  the  evaporative  cooling  which  results  maintains  the 
body  in  thermal  balance  with  the  environment.   If  the  body  demands  it, 
the  sweat  rate  will  increase  until  the  body  is  covered  with  moisture. 
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There  is  a  limit,  however,  to  how  much  moisture  the  air  will  hold, 
and  how  fast  it  will  take  it.   This  depends  upon  the  difference  in  the 
vapor  pressure  at  the  skin  surface  and  in  the  ambient  air  for  the  given 
body  surface  area.  As  the  ambient  temperature  decreases  and  the  rela- 
tive humidity  decreases,  the  air  will  accept  more  and  more  water-vapor. 

A  point  is  then  reached  where  the  air  can  take  no  more  moisture.   E 

r  max 

is  then  the  maximum  cooling  to  be  obtained  by  the  evaporation  of  sweat 
under  those  environmental  conditions.   A  higher  sweat  rate  will  pro- 
duce no  more  body  cooling  and  if  this  is  not  sufficient,  the  body  will 
store  heat.   As  can  be  seen  from  any  of  the  Figures  (JJ-17)  ,  if  the 
relative  humidity  can  be  maintained  at  a  low  value,  the  zone  of  evapor- 
ative regulation  is  considerably  wider  than  for  a  high  value  of  RH. 

The  effect  of  helium  on  the  comfort  zone  is  to  increase  its 
slope  and  make  it  narrower.   Increasing  pressure  has  the  same  effect 
since  equation  I  111  shows  that  an  increase  in  C  increases  the  influence 
of  each  degree  F  difference  between  skin  temperature  and  ambient  temper- 
ature. 

By  comparing  the  curves  for  one  atmosphere  helium  with  those  for 
helium  at  20  atmospheres,  it  can  be  seen  that  a  much  higher  ambient 
temperature  is  needed  at  twenty  atmospheres  to  offset  the  same  radiant 
temperature.   This  is  due  to  the  increase  in  convection  heat  transfer 
with  pressure  as  shown  in  Figure  (8)  .   While  convection  is  increasing, 
metabolism,  radiation  and  evaporation  remain  relatively  unchanged. 
It  is,  therefore,  necessary  to  bring  the  ambient  temperature  closer 
and  closer  to  skin  temperature  in  order  to  reduce  the  convective  heat 
loss.   Since  thin  walled  oceanographic  research  vessels  will  have 
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their  wall  temperatures  governed  by  the  temperature  of  the  water  which 

surrounds  them,  it  will  probably  become  necessary  in  the  larger  vessels 

to  construct  a  radiation  shield  between  the  walls  and  the  occupants. 

When  this  is  done,  the  T   .  will  approximate  T  . 

mrt       rr  a 

Comparing  the  Tft  -  Tmrt  line  on  the  graphs  for  one  and  twenty 
atmospheres  for  a  skin  temperature  of  92  F,  it  is  seen  that  this  line 
crosses  the  zone  of  minimal  evaporative  cooling  at  84  F  on  the  one 
atmosphere  chart  and  89  F  on  the  chart  for  twenty  atmospheres. 

While  a  one  atmosphere  100%  helium  environment  is  not  a  practical 
environment  in  which  to  live,  Table  IV  shows  that  at  a  dive  depth  of 
645  feet,  the  atmosphere  for  heat  transfer  purposes  is  100%  helium. 

Table  IV 
Composition  of  the  Atmosphere  vs  Depth  of  Dive 


Depth 

P 

ressure 

ft 

atm 

psia 

0 

1 

14.7 

32.3 

2 

29.4 

64.5 

3 

44.1 

96.8 

4 

58.8 

161.5 

5 

73.5 

290.5 

10 

147.0 

645 

20 

294.0 

°2 


He 


PO  psia  %  Vol   PRe  psia  %  Vol 

3.09  21     11.6  79 

3.09  10.5   26.3  89.5 

3.09  7.0   41.0  93 

3.09  5.25  55.7  94.75 

3.09  4.20  70.4  95.8 

3.09  2.1  143.9  97.9 

3.09  1.05  290.9  98.95 


Figure  3  shows  the  effect  of  this  dive  depth  on  the  modified  coefficient 
of  convection. 
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Figure  9.   Ambient  Temperature  vs  Mean  Radiant  Temperature 
Air  P  =  1  atm,  Tg  =  90F 

Test  Point  MRT  from  Page  51.  Numeral  indicates 
Test  Number. 
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Figure  10.   Ambient  Temperature  vs  Mean  Radiant  Temperature 
Air  P  =  1  atm,  Tg  =  92F 

Test  Point  MRT  from  Page  51.  Numeral  indicates 
Test  Number. 
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Figure  11.   Ambient  Temperature  vs  Mean  Radiant  Temperature 
Air  P  =  1  atm,  Tg  =  94F 

Test  Point  MRT  from  Page  51.  Numeral  indicates 
Test  Number. 
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Figure  12.   Ambient  Temperature  vs  Mean  Radiant  Temperature 


Air  P  =  1  atm.   Tg  =  96  F 
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Figure  13.  Ambient  Temperature  vs  Mean  Radiant  Temperature 
Helium  =  1  atm,  T  =  9 OF 

Test  Point  MRT  fr&n  Page  51.   Numeral  indicates 
Test  Number. 
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Figure  14.   Ambient  Temperature  vs  Mean  Radiant  Temperature 
Helium  =  1  atm,  Tg  =  92F . 

Test  Point  MRT  from  Page  51.   Numeral  indicates 
Test  Number. 
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Figure  16.   Ambient  Temperature  vs  Mean  Radiant  Temperature 
Helium  =  1  atm,  Tg  =  96F 

Test  Point  MRT  from  Page  51.   Numeral  indicates 
Test  Number. 
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Figure  17.   Ambient  Temperature  vs  Mean  Radiant  Temperature 


Helium  =  20  atm,  Tg  =  90-96F 
Test  Point  MRT  from  Page  51. 
Test  Number. 
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Experimental  Results 

Comfort  vote  data  are  shown  in  Table  V.   Investigators  experi- 
menting with  helium-oxygen  atmospheres  at  high  pressures  have  reported 
ambient  temperatures  adjusted  for  comfort  higher  than  in  an  environ- 
ment of  one  atmosphere  air  (30) (31) .   This  is  not  wholly  explained  by 
the  higher  conductivity  of  helium,  however,  since  Figure  (8)  shows  that 
the  convection  coefficient  for  air  also  increases  as  the  square  root 
of  the  pressure  ratio. 

Referring  to  Table  V,  in  two  tests,  practically  all  conditions 
were  identical  except  that  No.  3  was  in  helium  while  No.  8  was  in  air. 
Comparing  the  subjective  response,  it  is  seen  that  the  subject  was 
slightly  warm  in  air  while  he  was  comfortable  in  helium.   Although 
this  one  comparison  is  not  enough  to  form  a  conclusion,  it  is  an  indi- 
cation that  the  comfort  zone  in  helium  occurs  at  higher  values  of 
ambient  gas  temperatures. 

It  should  be  also  noted  that  the  mean  skin  temperature  in  test 
No.  3,  helium,  was  approximately  0.6  F  lower  than  in  No.  8  air.   This 
is  in  agreement  qualitatively  with  the  reports  of  Fox  et  al.  (4) 
and  Dianov  (32).   Tests  Nos.  1  and  7  appear  comparable,  yet  the  skin 
temperature  is  the  same  or  even  a  little  higher  in  helium.  It  should 
be  noted  that  the  subject  was  "cool"  in  both  cases  and  the  skin  was, 
therefore,  deprived  of  blood  due  to  constriction  of  blood  vessels  in 
an  attempt  to  conserve  body  heat.   Since  the  body  was  losing  too  much 
heat  in  either  case,  regardless  of  the  rate,  the  skin  temperature  was 
probably  lowered  an  equal  amount. 

The  experimental  conditions  corresponding  to  the  measured  skin 


Table  V 


Experimental  Comfort  Responses 
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Helium 


Test 

Number 

T„  (F) 

He 

RH 
% 

22 

MRT 
F 

75 

MRT„ 
F 

76.8 

T 
FS 

90.2 

Response 

1 

77.6 

Cool 

2 

82.5 

24 

75 

80.9 

91.7 

Comfortable 

3 

85.3 

33 

80 

82.8 

94.0 

Comfortable 

4 

88.2 

28 

75 

86 

94.6 

Warm 

5 

89.3 

34 

85 

86.8 

94.6 

Warm  -  Some  Sweating 

6 

92.6 

54 

88 

89.5 

95.9 

Warm  -  Profuse  Sweati 

where 


Air 


Test  Room 

7 

80 

30 

80 

89.9 

Cool 

Test  Cell 

8 
Test  Room 

85 

50 

80 

82.8 

94.6 

Slightly  warm 

9 

95 

30 

80 

94.8 

Warm 

TH   =  Helium  temperature ,  F 
RH   =  Relative  Humicity,  % 


T    =  Mean  skin  temperature,  F 


MRTj  =  Surface  temperature  of  test  room  walls,  F 

MRT2  =  Surface  temperature  of  polyethylene  test  cell,  F 
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temperatures  have  been  plotted  in  Figures  9  through  16  unadjusted 
for  the  difference  in  metabolic  rate  of  400  Btu/hr  and  the  subject's 
recorded  metabolic  rate.  Unadjusted,  the  position  of  the  test  points 
on  the  curves  correlate  very  well  with  the  subject's  comfort  response. 
There  is  reason  to  believe  that  the  subject's  metabolic  rate  was 
nearly  400  Btu/hr.   This  is  explained  later  when  the  heat  balance  is 
discussed.  The  results  of  the  heat  transfer  calculations  are  shown 
in  Table  VI.   The  experimental  analysis  was  done  for  the  most  part 
using  equations  stated  earlier  in  the  text. 

Radiation  heat  loss,  R,  was  evaluated  using  equation  |3|  p.  8 
and  was  also  evaluated  using  Gebhart's  method  (33)  for  enclosure 
problems.  Although  this  method  takes  into  account  radiation  which 
reaches  the  body  due  to  reflections  from  other  sources,  the  results 
were  not  significantly  different  from  those  obtained  by  using  equation  J  3 1 

The  convection  heat  loss  C  was  calculated  using  equation  |ll|  p.  16. 
The  total  DuBois  surface  area  was  used  for  A  .  The  experiment  described 
in  this  report  would  fall  into  the  free  convection  category  since  the 
subject  was  essentially  motionless  and  the  tube  leaking  helium  into  the 
test  cell  was  blocked  to  prevent  the  helium  from  blowing  on  the  sub- 
ject.  In  fact,  the  subject  was  unable  to  detect  the  flow  of  helium 
without  bringing  his  foot  within  inches  of  the  tube  opening.   The 
evaporative  heat  loss  was  calculated  by  the  following  formula. 

E  =  hfg  (Wf  -  W±)  -  (0.004355)  VC02  -  (0.003286)  VC02  Btu/hr   |l5| 
where 

h   =  heat  of  vaporization  at  mean  skin  temperature,  Btu/lb 

W^  =  subject  weight  after  1  hour,  lb 
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W    ■  subject  initial  weight,  lb 
i 

VCO  =  CCL  produced,  liters/hr 

V02  =  0  consumed,  liters/hr 

0.004355  =  weight  of  one  standard  liter  of  C02  in  lbs 

0.003286  =  weight  of  one  standard  liter  of  02  in  lbs 

For  a  respiratory  quotient  of  0.82,  the  equation  becomes: 

E  =  h   (w  -  W  )  -  (0.004355) (0.82)  -  (0.003286)  V0- 
fg   f    i  « 

=  hf  AW  -  (0.000285)  V02 
Also  since  the  subject's  mean  0„  consumption  was  15.2  1/hr,  the  equa- 
tion reduces  to: 

E  =  hf   AW  -  (0.000285) (16.2) 
=  hf   AW  -  0.0046,  Btu/hr 

The  measured  O2  consumption  was  corrected  to  standard  temperature 
and  pressure  (dry).   The  caloric  value  of  43  KCai/liter  is  equivalent 
to  0.0189  Btu/cc.   Since  the  02  consumption  was  measured  in  cc/min, 
this  was  multiplied  by  1.138  Btu/cc-hr  to  give  the  metabolic  rate  in  Btu/hr. 

In  the  discussion  which  follows,  the  metabolic  rate  will  be  repre- 
sented by  an  M,  the  radiative  heat  loss  by  an  R,  the  convective  heat 
loss  by  a  C  and  the  evaporative  heat  loss  by  an  E,  all  Btu/hr. 

When  M,  E,  R,  and  C,  were  calculated,  it  became  apparent  that  using 
the  existing  equations  the  metabolic  rate  would  not  equal  the  theoretical 
heat  loss,  although  the  storage  of  heat  was  assumed  to  be  zero.   In 
test  No.  6,  an  attempt  was  made  to  obtain  a  helium  temperature  equal  to 
the  skin  temperature.  The  test  room  walls  were  maintained  at  95  F  and 
the  floor  at  75  F.  According  to  equation  |3 |  ,  this  would  give  a  value 
for  R,  =  85  Btu/hr  which  when  subtracted  from  a  metabolic  rate  of  325 
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would  leave  240  for  E  which  the  body  would  be  able  to  handle.   At 

Ta  =  93  F,  the  sweat  rate  was  such  that  a  further  increase  in  Ta  would 

have  caused  the  sweat  to  run  off  so  T  was  lowered  0.5  F.   When  the  data 

a 

were  analyzed,  it  was  found  that  both  the  calculated  C  and  R  exceeded  M-E. 

An  experiment  was  run  with  air  in  the  test  cell,  No.  8,  with  the 
conditions  of  No.  3  in  hopes  that  the  change  from  helium  to  air  would 
not  only  affect  the  comfort  vote,  but  would  bring  the  calculated  heat  ex- 
change of  the  subject  with  his  environment  to  zero. 

Two  tests,  Nos.  7  and  9,  were  then  run  with  the  same  subject  in 
the  same  equipment  only  outside  the  test  cell.   The  heat  balances  then 
checked. 

By  achieveing  T  =  T  in  Test  No.  9,  the  convection  term  was  zero 
a    s 

and  M-E-R  =  0. 

The  evaporative  weight  loss  was  measured  on  a  scale  graduated  in 
increments  of  0.01  lbs.  Therefore,  an  accuracy  of  ±0.01  lbs  is  rea- 
sonable.  This,  when  multiplied  by  an  average  heat  of  vaporization  of 
10A0  Btu/lbw  equals  an  error  of  ±10.4  Btu.   The  evaporative  heat  losses 
versus  skin  temperatures  for  the  tests  are  plotted  in  Figure  18  with 
similar  data  reported  by  Gagge  (16).   The  agreement  is  good  and  lends 
credence  to  the  experimental  measurements  of  weight  loss.   One  expla- 
nation for  the  two  low  readings  is  that  the  short  time  period  of  the 
experiment  did  not  allow  the  "corneum  stratum",  horny  layer  of  the 
skin,  to  become  saturated,  and  diffusion  of  water  vapor  was  inhibited, 
Fanger  (34) .   This  also  shows  that  the  evaporative  weight  loss  is  inde- 
pendent of  the  gas  to  which  the  skin  is  exposed  and  is  a  function  of 
skin  temperature. 
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It  may  also  be  noted  (Figure  18)  that  as  the  need  for  more  evapora- 
tive cooling  arises,  the  skin  temperature  also  increases  due  to  vaso- 
motor regulation.  Thus,  in  terms  of  the  comfort  lines  in  Figures  8 
through  16  seldom,  if  ever,  will  F^^  be  reached  for  Tg  =  90  or  92  F. 
When  the  environmental  conditions  are  such  that  larger  amounts  of 
evaporative  cooling  are  required  the  skin  temperature  will  have  risen 
to  94  or  96  F. 

The  measured  metabolic  rate  appears  low  when  compared  with  values 
from  the  literature,  (5) (6) .  However,  the  subject's  metabolic  rate 
while  seated  and  sedentary  was  measured  seven  times  while  in  the  test 
cell  and  at  least  five  times  outside  the  test  cell.   The  equipment  had 
been  used  previously  by  the  same  operator  with  good  results  and  was 
considered  in  excellent  condition.   As  can  be  seen  in  Table  VI  the 
differences  among  the  tests  are  no  greater  than  would  be  expected. 
The  subject's  basal  metabolic  rate  measured  at  the  Kansas  State  Univer- 
sity Student  Health  Center  was  found  to  be  280  Btu/hr  which  is  within 
the  normal  range.   It  was  suggested  that  since  the  subject  completely 
relaxed  during  the  time  the  metabolic  rate  was  being  measured,  this 
rate  was  not  the  same  as  it  would  have  been  had  he  continued  to  measure 
his  weight,  operate  the  relative  humidity  sensor  and  record  his  pulse 
and  subjective  reactions  as  he  did  otherwise.  A  metabolic  rate  was 
measured  for  the  subject  while  typing  and  was  found  to  be  400  Btu/hr. 
This  is  .over  than  the  550  Btu/hr  suggested  for  this  activity  level  (5). 
However,  no  typing  rate  was  recorded  in  either  case.   If  the  subject 
was  assumed  to  have  had  a  metabolic  rate  of  400  Btu/hr,  for  45  minutes 
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and  a  metabolic  rate  of  325  Btu/hr  for  15  minutes,  his  combined  meta- 
bolic rate  would  have  been  381  Btu/hr.   This  does  not  bring  the  heat 
exchange  into  balance,  but  it  would  influence  the  correction  of  the 
experimental  points  plotted  on  Figures  9  through  16.   Since  this  dis- 
crepancy existed  between  expected  and  measured  metabolic  rates,  the 
calculations  were  checked  for  errors.  None  were  found  and  the  dis- 
crepancy remains  unexplained. 

In  test  No.  9,  M-E-R  =0,  and  there  is  substantial  reason  to  believe 
M  and  E  are  valid.   Therefore,  the  value  for  R  and  its  method  of  cal- 
culation should  also  be  valid.   In  test  No.  8  a  convection  term  is 
added  and  the  energy  produced  by  the  subject  equals  the  energy  lost. 

Once  the  subject  enters  the  test  cell,  however,  the  equations  no 
longer  check.  The  heat  loss  due  to  radiation  exchange  is  now  compli- 
cated by  the  presence  of  the  polyethylene  test  cell.   As  mentioned 
earlier,  this  material  when  clean  and  dry  will  transmit  85%  of  the  infra- 
red radiation  incident  upon  it.   Therefore,  one  solution  was  to  consider 
the  test  cell  transparent  and  keep  the  shape  factor  F  an_fioor  =  0.33. 
This  is  MRT  in  Table  V  and  results  in  a  radiation  heat  exchange  R-,  . 
The  second  approach  was  to  consider  the  test  cell  covered  on  the  inside 
with  a  thin  film  of  moisture  which  it  had  acquired  while  in  the  con- 
figuration shown  in  Figure  3,  collapsed.  The  shape  factor  Fman_f  -,,  Qr~  0-15 
and  the  wall  temperature  of  the  test  cell  replaces  the  wall  temperature 
of  the  test  room  in  the  calculation  of  MRT„   This  resulted  in  MRT2  and  R2 
in  Table  V.   As  can  be  seen  by  comparing  the  M-E-R„  column  with  the 
theoretical  convection  column,  the  heat  transfer  equation  did  not  balance. 

A  complicated  radiation  heat  exchange  analysis  was  attempted 
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assuming  a  transmissivity  of  (0.85)  for  infrared  radiation  incident 

upon  the  test  cell.   This  resulted  in  a  higher  radiation  heat  loss  and 
a  greater  heat  balance  discrepancy  than  did  the  other  solutions.   The 
only  factor  in  the  radiation  term  which  could  change  besides  the  MRT 
would  be  the  radiation  area  of  the  body  and  this  seems  fixed. 

The  convection  heat  transfer  equation  is  affected  by  T  ,  T  ,  A  , 
a  characteristic  dimension  and  the  thermal  properties  of  the  ambient 
gas.  Mean  skin  temperatures  were  in  the  normal  range  of  90-96  F. 
Due  to  the  fact  a  Tg  calculated  from  equation  13  will  vary  with  indi- 
vidual shapes,  any  error  in  the  measured  and  calculated  skin  temper- 
ature is  probably  no  greater  than  the  error  inherent  in  the  assumptions 
made  to  obtain  the  value.  Since  the  subject  is  of  average  height  and 
weight  this  error  should  be  small.  By  placing  the  thermisters  at 
approximately  the  same  location  for  each  test  and  monitoring  them  every 
ten  minutes  of  each  test,  it  was  quickly  noticed  when  a  thermister 
was  loose  and  exposed  to  the  atmosphere  or  was  being  pressed  between 
the  skin  and  the  chair  and  thereby  was  measuring  an  invalid  skin  tem- 
perature. 

The  surface  area  of  the  subject  is  a  factor  also  to  be  considered. 
DuBois  (8)  claims  his  height-weight  formula  is  accurate  to  one  decimal 
place  or  has  an  error  less  than  2%.   In  the  sitting  position,  the 
convection  area  should  be  within  the  2%  error  of  the  total  body  sur- 
face area. 

The  lower  skin  temperature  in  the  one  helium  test  compared  to  its 
air  counterpart  shows  that  the  higher  convection  loss  in  helium  does 
exist.   At  this  point  there  does  not  appear  to  be  a  sound  reason  why  the 
equation  j 13 |  should  not  give  valid  results. 
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RECOMMENDATIONS 

In  any  project  of  this  nature  a  lot  of  time  is  taken  finding  out 
the  who,  what,  where,  and  how  of  assembling  equipment  and  conducting  the 
research.  Hindsight  is  greater  than  or  equal  to  foresight,  and  the 
following  suggestions  could  be  considered  part  of  the  "delta  experience" 
gained  from  this  work. 

To  avoid  the  problems  with  radiation  heat  transfer  reported  in 
this  paper,  two  approaches  are  suggested.   (1)   Coat  the  test  cell  with 
a  flexible  paint.  This  would  assure  the  experimenter  that  the  sub- 
ject was  radiating  to  the  test  cell  alone.   (2)  Purchase  a  custom- 
made  test  cell  with  known  thermal  radiation  properties. 

Silica  gel  could  be  used  inside  the  test  cell  to  control  the  rela- 
tive humidity.   The  floor  temperature  could  be  maintained  equal  to  the 
test  cell  wall  temperature.   This  would  minimize  the  radiation  shape 
factor  difficulties. 

Although  the  evaporative  weight  losses  appear  to  be  reasonable, 
a  better  method  of  measurement  should  be  devised.   It  would  be  better 
to  have  someone  other  than  the  subject  record  the  weights.   Although 
if  there  is  any  chance  the  subject's  position  may  affect  the  reading, 
then  the  subject  must  be  instructed  to  find  a  reference  position  and 
assume  this  position  whenever  a  weight  is  taken.   The  best  system  would 
be  to  continuously  record  the  body  weight.   This  could  be  done  with  a 
device  employing  strain  gauges  and  a  strip  chart  recorder. 

The  mask  and  valve  are  cumbersome  and  somewhat  uncomfortable.   How- 
ever, they  are  far  superior  to  a  mouthpiece  and  noseclip  combination, 
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especially  for  a  test  of  one  hour  or  longer.   Both  the  mask  and  mouth- 
piece cannot  be  eliminated  without  going  to  a  mixture  of  helium  and 
oxygen . 

More  subjects  would  be  desirable  to  allow  a  statistical  analysis 
of  the  results.   The  subjects  should  be  thoroughly  calibrated  in  air 
for  metabolic  rates  at  various  activities,  mean  skin  and  rectal  temper- 
atures, sweat  rates  and  comfort  zones.   Several  trials  in  air  at  various 

conditions  includint  T  =  T  and  T  =  T    should  be  made  and  heat 

s    a      s    mrt 

balances  calculated. 

After  this  calibration  any  irregularities  which  occur  in  the 

actual  test  should  be  quickly  noticed. 

Two  sets  of  test  under  the  conditions  T  =  T  and  T  =  T 

s    a      s    mrt 

should  be  conducted  and  carefully  analyzed  to  try  to  determine  coeffi- 
cients of  radiation  and  convection  heat  transfer  before  any  other  tests 
are  attempted. 

Several  tests  at  each  condition  with  each  subject  should  be  made. 
The  subject  should  be  exposed  to  the  test  conditions  one  and  one-half 
hours  before  a  steady  state  experiment  is  conducted,  (in  the  comfort 
zone) . 

The  mass  spectrometer  in  the  Biochemistry  department  can  be  used 
to  obtain  a  gas  analysis  which  will  show  the  percentage  of  helium  and 
the  percentage  of  air  in  a  sample.  The  only  person  qualified  to  per- 
form this  analysis  left  Kansas  State  University  the  week-end  preceding 
the  first  test  in  helium.   The  biochemistry  personnel  are  very  helpful 
and  a  gas  analysis  should  be  made  of  each  test. 
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MONS$      JOB   TMRT  HELIUM  OR  AIR   W/  -E+/-C+/-R 
MON$$      COMT  20,10.     G.  SMITH     MECH.  GEPI. 
MO.\*i     ASGN  WJB.12 
MOMS!      ASG.W  mG0,1A 

t$$     MODE  GU.TEST 
MOMS*      EXEQ  FORTRAN,  ,,,,,,  SOLVE 

REAL  H 

DIMENSION  TS(6) , TA( 11) ,WU{2) ,RH(4) ,A(6) 

CO MMON  C I , RL , PR , AB, Gt L 

1  FORMAT  (4F18.6) 

2  F0RMA1  (SriO-l) 

3  FU  (MAT  (2FlQ.lt  I  LO) 

5  FORMAT  (4FL0.  1  ,«8i>FlO.  I  )  ,.., 

6  FORMAT  (  1H1,40X,2;:M=,F4.0,6X,2HP=,F8.3,6X,  6HHtLIUM///) 

7  FORMAT  (3lX,2HTi,,i0X,2HTA,9X,2HWU,9X,2HRH,  ;X,1M£,  11X.1HC, 
8DX,4HTMRT) 

0  FORMAT  (22X,2F12.2,F11.?,F10.2,2F12.2,F13.2) 

RFADI 1,3)M,P,L 

READ(1,2)TS 

READ  ( 1,2)  I  A 

READ  { 1,2) WU 

READ  ( 1,2)KH 

READl  1,  DA 

READ  { i,5)  C1,RL,PR,AB, 

READ  ( 1,2)  SIG,AR,EPS,C. 

DO  41=1,4 

WRITE  <3,6)M,P 

WRITE  (3,7) 

TS  I  =  TSID+460. 

TSI2  =TSI*TSI 

TSJ  =  (TSI  I)+459.69)/l0 •'■. 

VPS=A( U+AI2)* IS1+A{3)*!SI2+A(4)*TSI*TSI2+A(5)*TSI2*TSI2+A(6)» 
BTSI*TSI2»TSI2 

DU  413=1,11 

CALL  C  CALC  (TSU ),TA(  I  3 ) , TF , P ,CC, C ) 

TAI=TA{ I3)+46u. 

TAI2=TAI»TAI 

VPA=A{ D+A(2)*TAI+A(3)*TAI2+A(4)«TAI*TAl2+A{5)*TAI2»rAl2+A(6)* 
BTAI*TAI2*TAI2 

IFITAI I3).LE.TS{I) )  GO  TO  10 

c=-c 

10  DO  411=1,2 
DO  412=1,4 

IF(VPS-RH{ I2)»VPA)11,11,12 

11  E=.  . 

GU  TO  13 

12  E  =  WU( I l)*AQ»CF»(VPS-RH(  I2)*VPA) 

13  TMRT  =(TSJ»*4- (H-C-E)/(-SIG*AR*EPS)  )**.25 
TMRT  =(TMRT*10( ,. ) -459. 69 

4  WRI TE{3,B)TS( I ),TA( 13) ,WU(I1) ,RH(I2) ,E,C,TMRT 
STOP 

END 
MON$*      EXEQ  FOR  IRAN 
SOUROOTIMECCALCITSl, TE1,TF,P,CC,C) 
C0MM0NC1 i RL , PR , AB , G , L 
TS=TSl+459.69 
TE=T£l+459.69 
TF=.5MTS+TE) 
BE1 A=1./TF 
GO  fO(7,8),L 

7  RH0=P/(2.6  829*TF+P*RETA*».25*(  .3059-1  .  845»f,QRT  (  BETA  )-.  fi  22  *rtET  A  )  ) 
RMU=.8315*.1E-02»TF**.647 

RK=(7920.+9.62*(  TF-459  .  69  )  +  9.  *  {  P/14.  696-1 .  )  )*.lE-04 
GO    TO    10 

8  PR=HR+.031^ 
RHO=P/{0.3  70484«TF) 
TFK=0.555*TF 

RMU=( (TFK*SCRT( TFK) )/( TFK+110.4) )*3.527E-0  3 

R!<=(?42.*.6325*SQRT(TFK)*I.E-05)/{  1 . +245.4*10 .»*( -12 ./ TFO / TFK ) 
10    CC=C1*RK*(G*PR*BEIA*RHU*RH0/(RL*RMU*RMU) )**.25 
C=AB*CC*(ABS(  TS-Tfc )  )**1.25 
TF= TF-459. 69 

5  RETURN 
END 
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Explanation  of  Fortran  Program 

This  program  is  written  for  the  IBM  1410  computer  using  the  for- 
tran  language  found  in  the  IBM  manual  File  No.  1410/7010-25  Form 
C28-0328-2. 

Before  attempting  to  use  this  program  enough  of  the  theoretical 
analysis  presented  in  the  text  of  this  report  should  be  read  to  familiar- 
ize the  user  with  the  nomenclature. 

The  program  solves  the  heat  balance  equation,  No.  13  in  the  text 
of  this  report,  which  reads,  metabolic  rate  (M)  =  evaporative  heat 
loss  (E)  ±  convective  heat  loss  C  ±  radiative  heat  loss  R.  The  meta- 
bolic rate  M  is  read  into  the  machine.   The  evaporative  heat  loss,  E, 
is  calculated  in  statement  No.  12.   The  convective  heat  loss,  C,  is 
calculated  in  the  subroutine  named  C  CALC.   The  radiation  term  has 
not  been  calculated  separately.   Instead,  equation  13  has  been  solved 

for  T    as  is  shown  in  statement  13  of  the  program. 

mrt  r  ° 

If  the  program  is  used  on  the  IBM  1410  using  the  PR-155  system, 
the  only  cards  needed  are  the  data  cards . 

The  first  data  contains  M,  metabolic  rate  in  Btu/hr,  P, 
pressure  in  psia  and  L,  a  floating  point  number  which  when  a  1 
designates  helium  and  a  2  designates  air.   The  next  data  card  contains 
the  skin  temperature  (Tg)  in  degrees  F.  The  third  data  card  contains 

the  ambient  temperatures  in  degrees  F.   The  fourth  data  card  lists 

2 
the  values  for  Gagge's  wetted  area  in  KCal/cm  hr-m  and  the  fifth 

values  for  relative  humidity.   The  constants  A  listed  on  the  sixth  data 

card  are  used  to  generate  values  of  vapor  pressure  in  psia  from  values 

of  temperature  in  degrees  F:   A(l)  =  -0.7321633  •  10~5 ,  A(2)  = 


69 

_5 
0.0604572,  A(3)  =  0.0021633962,  A(4)  =  -0.7321633  -10  ,  A(5)  = 

0.84661484  '10"  .   On  the  seventh  data  card  are  contained  values 

for:   the  constant  listed  as  CI,  dimensionless,  the  significant 

dimension  L,ft,  the  Prandtl  number,  dimensionless,  the  total  body 

9 
surface  area,  square  feet  and  the  acceleration  of  gravity  in  ft/hr  . 

In  the  text  the  values  for  CI  and  L  are  combined  in  one  number  but 

when  the  program  was  written  a  value  of  CI  was  selected  as  0.555  and 

0  25 
the  value  of  L  altered  to  fit  selected  values  for  <j>  =  Cl/L  '   .   The 

final  data  card  contains  the  values  of  the  Stefan  Boltzman  Constant, 

the  radiation  area  and  emissivity  of  the  skin,  and  the  conversion 

factor  1.90633  Btu-cm  Hg-m2/KCal-psia-f t2. 
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ABSTRACT 

Six  tests  were  conducted  involving  a  seated  nude  subject  exposed 
to  a  helium  environment  using  the  equipment  and  techniques  described 
in  the  body  of  this  report.   In  addition,  a  test  was  run  in  the  test 
cell  used  to  contain  helium  with  air  as  the  ambient  gas  and  two  tests 
were  run  in  air  using  the  same  equipment  with  the  exception  of  the 
polyethylene  bag.   During  these  tests  the  subject  recorded  his  sub- 
jective thermal  comfort.   Also  the  parameters  necessary  to  determine 
the  subject's  heat  balance  were  measured  and  calculations  were  made 
according  to  heat  transfer  theory. 

For  the  conditions  of  this  :.est,  the  subject  was  comfortable  with 
ambient  helium  temperatures  of  82-86  F. 

The  mean  skin  temperature  in  a  helium  atmosphere  at  85  F  was  0.6  F 
lower  than  the  skin  temperature  in  85  F  air. 

The  evaporative  heat  loss  was  found  to  vary  with  the  skin 
temperature  and  is  independent  of  the  ambient  gas  for  a  given  skin 
temperature. 

A  theoretical  analysis  was  made  for  the  heat  transfer  from  the 
human  body  according  to  heat  transfer  theory.   The  theory  determines 
a  mathematical  equation  for  heat  transfer  and  heat  balance  in  the 
zone  of  thermal  equilibrium  and  comfort.   The  results  of  the  equation 
are  plotted  on  a  graph  of  mean  radiant  temperature  versus  ambient 
temperature  at  pressures  of  one  and  twenty  atmospheres  and  experimen- 
tal points  are  compared  with  theoretical  predictions  for  one  atmosphere 
pressure. 
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